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The photon spectrum in B — ♦ X s "/ decay, where X s is any strange hadronic state, is studied 
using a data sample of 88.5 x 10 6 e + e~ — > X(4S) — > BB decays collected by the BABAR experiment 
at SLAC. The partial branching fraction, AB(B -> X s j) = (3.67 ± 0.29(stat.) ± 0.34(sys.) ± 
0.29{model)) x 10~ 4 , the first moment (E 7 ) = 2.288 ± 0.025 ± 0.017 ± 0.015 GeV and the second 
moment (E%) = 0.0328 ± 0.0040 ± 0.0023 ± 0.0036 GeV 2 are measured for the photon energy range 
1.9 GeV < E 1 < 2.7 GeV. They are also measured for narrower E-, ranges. The moments are then 
fit to recent theoretical calculations to extract the Heavy Quark Expansion parameters, mj, and 
H%, and to extrapolate the partial branching fraction to Ej > 1.6 GeV. In addition, the direct CP 
asymmetry Acp(B —* X s+ d^) is measured to be —0.110 ± 0A15(stat.) ± 0.017(sys.). 



PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 

In the Standard Model (SM) the radiative decay of 
the b quark, b — > s~/, proceeds via a loop diagram, 
and is sensitive to possible new physics, with new heavy 
particles participating in the loop 0. Next-to- leading- 
order SM calculations for the branching fraction give 
B(B -> X sl ) = (3.6ll°:f£) x 10- 4 (£; 7 > 1.6 GeV) 0, 
and calculations to higher order, which are expected to 
considerably decrease the uncertainty, are currently un- 
derway 0. The shape of the photon energy spectrum, 
which is insensitive to non-SM physics [4j, can be used to 
determine the Heavy Quark Expansion (HQE) parame- 
ters, m\y and fi^ [5j,|fJ, related to the mass and momentum 
of the b quark within the B meson. These parameters 
can be used to reduce the error in the extraction of the 
CKM matrix elements \V c b\ and \V u b\ from semi-leptonic 
B-meson decays 0. New physics can also significantly 
enhance the direct CP asymmetry for b —> sj and b — » d^j 
decay 0, A CP = ^^+^7)-r(^7+^7) which is 

Ri 10~ 9 in the SM [8j. Measurements of this joint asym- 
metry complement those of Acp in b — > 57 9] to con- 
strain new physics models. 

This letter reports on a fully-inclusive analysis of B — > 
A" s 7 decays collected from e + e~ — > T(4S) — ► BB, where 
the photon from the decay of one B meson is measured, 
but the X s is not reconstructed. This avoids incurring 
large uncertainties from the modeling of the X s fragmen- 
tation, but at the cost of high backgrounds which need to 
be strongly suppressed. The principal backgrounds are 
from other BB decays containing a high energy photon 
and from continuum qq (q = udsc) and t + t~ events. The 
continuum background, including a contribution from 
initial state radiation (ISR), is suppressed principally by 
requiring a high-momentum lepton from the non-signal B 
decay, and also by discriminating against its more jet-like 
topology. The BB background to high energy photons, 
dominated by ir° and 77 decays, is reduced by vetoing 
on reconstructed ir° or r\. The residual continuum back- 
ground is subtracted using off-resonance data taken at a 
center-of-mass energy 40MeV below that of the T(4S), 
while the remaining BB background is estimated using 
a Monte Carlo simulation which has been checked and 



corrected using data control samples. Previous inclusive 
measurements of B — > X s j have been presented by the 
CLEO [13, BELLE 11] and BABAR collaborations 
using alternative techniques which incur different system- 
atic uncertainties. 

The results presented are based on data collected 
with the BABAR detector [l3| at the PEP-II asymmetric- 
energy e + e~ collider located at the Stanford Linear Ac- 
celerator Center. The on-resonance integrated lumi- 
nosity is 81.5 fb~ , corresponding to 88.5 million BB 
events. Additionally, 9.6 fb -1 of off-resonance data are 
used in the continuum background subtraction. The 
BABAR Monte Carlo simulation program, based on 
GEANT4 Q, EVTGEN 13 and JETSET 0, is used 
to generate samples of B + B~ and B°B° (excluding sig- 
nal channels), qq, t + t~, and signal events. The signal 
models used to calculate efficiencies are based on refer- 
ences ("kinetic scheme") and 0] ("shape function 
scheme" ) and on an earlier calculation [4j ( "KN" ) . These 
predictions approximate the X s resonance structure with 
a smooth distribution in mx s ■ This is reasonable except 
at the lowest masses where the if* (892) dominates the 
spectrum. Hence the portion of the rnx s spectrum be- 
low l.lGeV/c 2 is replaced by a Breit-Wigner K*(&92) 
distribution. The analysis was done "blind" in the range 
of reconstructed photon energy E* from 1.9 to 2.9 GeV 
(asterisk denotes the T(4S) rest frame); that is, the on- 
resonance data were not looked at until all selection re- 
quirements were set and the corrected backgrounds de- 
termined. The signal range is limited by high BB back- 
grounds at low E*. 

The event selection begins by finding at least one pho- 
ton candidate with, 1.6 < E* < 3.4 GeV, in the event. A 
photon candidate is a localized electromagnetic calorime- 
ter energy deposit with a lateral profile consistent with 
that of a single photon. It is required to be isolated by 
25 cm from any other energy deposit and to be well con- 
tained in the calorimeter (—0.74 < cos6* 7 < 0.93), where 
9~f is the polar angle with respect to the beam-axis. Pho- 
tons that are consistent with originating from an identifi- 
able 7T° or i] — > 77 decay are vetoed. Hadronic events are 
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selected by requiring at least three reconstructed charged 
particles and the normalized second Fox- Wolfram mo- 
ment i?2 to be less than 0.55. To reduce radiative Bhabha 
and two-photon backgrounds, the number of charged par- 
ticles plus half the number of photons with energy above 
0.08 GeV is required to be > 4.5. 

Event shape variables are used to exploit the difference 
in topology of isotropic BB events and jet-like continuum 
events. This is accomplished by the R 2 requirement as 
well as a single linear discriminant formed from nineteen 
different variables. Eighteen of the quantities are the sum 
of charged and neutral energy found in 10-degree cones 
(from to 180 degrees) centered on the photon candi- 
date direction; the photon energy is not included. Ad- 
ditionally the discriminant includes i? 2 /i?2, where R' 2 is 
the normalized second Fox- Wolfram moment calculated 
in the frame recoiling against the photon, which for ISR 
events is the qq rest frame. The discriminant coefficients 
were determined by maximizing the separation power be- 
tween simulated signal and continuum events. 

Lepton tagging further reduces the backgrounds from 
continuum events. About 20% of B mesons decay semi- 
leptonically to either e or fi. Leptons from hadron de- 
cays in continuum events tend to be at lower momentum. 
Since the tag lepton comes from the recoiling B meson, 
it does not compromise the inclusiveness of the B — > X s ^ 
selection. The tag lepton is required to have momentum 
p* e > 1.25GeV/c for electrons and p* > 1.5GeV/c for 
muons. Additionally requiring the photon-lcpton angle, 
cos6**^ > —0.7 removes more continuum background, in 
which the lepton and photon candidates tend to be back- 
to-back. Finally the presence of a relatively high-energy 
neutrino in semi-leptonic B decays is exploited by requir- 
ing the missing energy of the event, E^ iss > 0.8GeV/c. 
Virtually all of the tagging leptons arise from the decay 
B — > XJtv. The rate of such events in the simulation is 
corrected as a function of lepton momentum 0. 

The event selection is chosen to maximize the statis- 
tical significance of the expected signal using simulated 
signal (KN with m b =4.80 GeV/c 2 , = 0.30 GeV 2 ) and 
background events, allowing for the low statistics of the 
off-resonance data used for the subtraction of continuum 
background. After selection the low energy range, 1.6 < 
E* < 1.9 GeV, is dominated by the BB background, 
while the high energy range, 2.9 < E* < 3.4 GeV, is 
dominated by the continuum background; they provide 
control regions for the BB subtraction and continuum 
subtraction, respectively. The signal region lies between 
1.9 GeV and 2.7 GeV. The signal efficiency (» 1.6% for 
this E* range) depends on E* and the signal model, but 
has negligible dependence on the details of the fragmen- 
tation of the X s . 

The BB background is estimated with the simulated 
BB data set. It consists predominantly of photons orig- 
inating from 7T° or 77 decays (« 80%). Other significant 
sources are n's which fake photons by annihilating in the 



calorimeter and electrons that are misreconstructed or 
lost, or that undergo hard Bremsstrahlung. The 7t°(t/) 
background simulation is compared to data by using the 
same selection criteria as for B — > X s j but removing the 
Tr°(r]) vetos. The photon energy and lepton momentum 
thresholds are relaxed to E* > 1.0 GeV, p* > l.OGeV/c, 
jj* > 1.1 GeV/c to gain statistics. The yields of 7T°(?7) are 
measured in bins of E* ^ by fitting the 77 mass dis- 
tributions in on-resonance data, off-resonance data and 
simulated BB background. Correction factors to the tt° 
(r)) components of the BB simulation are derived from 
these yields, including a small adjustment for the dif- 
ferent efficiencies of the ir° (77) vetoes between data and 
simulation. As no n control sample could be isolated, this 
source of BB background is corrected by comparing in 
data and simulation the inclusive p yields in B decay and 
the calorimeter response to p's, using a A — > pn + sam- 
ple. The electron component of the BB simulation is cor- 
rected with electrons from a Bhabha data sample, taking 
into account the lower track multiplicity of these events 
compared to the signal events. Finally, the small contri- 
butions from u> and rf decays are corrected using inclusive 
B decay data. After including all corrections and system- 
atic errors the expected background yield from the simu- 
lation in the BB control region (1.6 < E* < 1.9 GeV) is 
1667 ± 54 events, compared to 1790 ± 64 events observed 
in data after continuum subtraction. Note that a small 
contribution in this region from the expected signal (« 20 
to 40 events) has been neglected in this comparison. In 
the high energy control region 2.9 < E* < 3.4 GeV the 
expected background is 390 ± 20 events, compared to 
393 ± 58 events observed in data. 




2.2 2.4 2.6 2.8 
E* y (GeV) 

FIG. 1: The photon energy spectrum after background sub- 
traction, uncorrected for efficiency. The inner error bars are 
statistical and the outer include systematic errors added in 
quadrature. The histograms show the spectra for values of 
mt and ^ from the best fits to the moments in the kinetic 
scheme (dashed) and shape function scheme (dotted), nor- 
malized to the data in the signal region. 

Figure ^ shows the measured spectrum for signal and 
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TABLE I: The measured partial branching fraction, first and second moment (±stat. ± syst. ± model) for different ranges of 
E-y in the B rest frame. 



E 1 (GeV) 


AB{B -► X sl ) (10- 4 ) 


(E,) (GeV) 


(E^) - (E,r (GeV^) 


1.9 to 2.7 

2.0 to 2.7 

2.1 to 2.7 

2.2 to 2.7 


3.67 ±0.29 ±0.34 ±0.29 

3.41 ±0.27 ±0.29 ±0.23 
2.97 ±0.24 ±0.25 ±0.17 

2. 42 ±0.21 ±0.20 ±0.13 


2.288 ± 0.025 ± 0.017 ± 0.015 
2.316 ± 0.016 ± 0.010 ± 0.013 
2.355 ± 0.014 ± 0.007 ± 0.011 
2.407 ± 0.012 ± 0.005 ± 0.008 


0.0328 ± 0.0040 ± 0.0023 ± 0.0036 
0.0266 ± 0.0026 ± 0.0010 ± 0.0020 
0.0191 ± 0.0019 ± 0.0006 ± 0.0015 
0.0116 ± 0.0014 ± 0.0004 ± 0.0005 



control regions after the BB and continuum backgrounds 
have been subtracted. To extract partial branching frac- 
tions (PBFs) and first and second moments from this 
spectrum it is necessary to first correct for efficiency. 
Theoretical predictions are made for the true E 1 in the 
B meson rest frame, whereas the experimental measure- 
ments are made with reconstructed E* in the T(AS) 
frame. Hence it is also necessary to correct for smear- 
ing due to the asymmetric calorimeter resolution and 
the Doppler shift between the T(4S) frame and the B 
rest frame. The efficiency and smearing corrections de- 
pend upon the assumed signal model (underlying the- 
ory and parameter values). In a broad selection of sig- 
nal models it is found that the efficiency for each E* 
range has a model-independent linear relationship to the 
mean E* in that range. Hence a nominal signal model 
is chosen for which the mean matches the data, and a 
model-dependence uncertainty is assigned to the PBFs 
and moments based on signal models within one (sta- 
tistical and systematic) standard deviation of the mea- 
sured mean E*. To correct for resolution smearing a 
small multiplicative correction to the PBF and small ad- 
ditive corrections to the first and second moments are 
computed using the nominal signal model, and an un- 
certainty assigned based on a conservative range of mod- 
els. The model-dependence uncertainty from the smear- 
ing correction is fully correlated with the corresponding 
uncertainty of the efficiency correction. 

The results for four energy ranges are given in Ta- 
ble 1 along with the statistical, systematic and model 
errors. The PBFs have been corrected to exclude a 
(4.0 ± 0.4)% H03 contribution from b -> dry. The sys- 
tematic errors are described below and the associated 
correlation matrices are given in the appendix. 

The most significant systematic uncertainty in the 
measurement of the spectrum is from the uncertainty 
in the corrections to the BB background simulation. It 
is due mostly to the statistical uncertainty on the cor- 
rection factors derived from the ir°(r)) control sample. 
The BB corrections depend on E*; the resulting corre- 
lations between the 100 MeV E* bins have been taken 
into account in the computation of the total systematic 
uncertainty in the PBFs and moments. For example, 
for 2.0 GeV < E 1 < 2.7 GeV, the BB corrections con- 
tribute 5.5% to a total systematic uncertainty of 8.5% of 
the PBF, and 0.008 GeV and 0.0009 GeV 2 of the total 



systematic uncertainty of the first and second moments, 
respectively. Additional contributions to the PBF un- 
certainty (added in quadrature), all energy-independent, 
come from the photon selection (3.3%) due to the pho- 
ton efficiency, determined with 7r°'s from r decay, and the 
isolation requirement, calorimeter energy scale and reso- 
lution, determined from B — > K*j decays and photons 
from virtual Compton scattering; efficiency of the event 
shape variable selection (3%), determined from a tt° con- 
trol sample; the semi-leptonic corrections (3%); lepton 
identification (2%) and the modeling of the X s fragmen- 
tation (1.5%). Additional uncertainties to the first and 
second moment, added in quadrature, come from the un- 
certainty in the calorimeter energy scale (0.006 GeV) and 
resolution (0.0004 GeV 2 ), respectively. 

The parameters mi, and /^ 2 , which are defined differ- 
ently in the kinetic (K) and shape function (SF) schemes, 
can be extracted by fitting theoretical predictions to the 
measured moments. The first moments for E 7 > 1.9 and 
2.0 GeV and the second moment for E 1 > 2.0 GeV are 
fitted, taking into account the correlations between the 
measured moments. As the moments are dependent on 
the assumed signal model due to the efficiency and res- 
olution smearing corrections, the signal model and the 
model-dependence errors are adjusted based on the re- 
sults of the fit and the moments are recomputed and refit. 
Only a few iterations are required until the result is sta- 
ble. In the kinetic scheme m b(K) = ^A^^-oaa GeV/c 2 
and nl {K) = 0.64i°;^t°; 2 ^GeV 2 , with a correlation of 
—0.93. The first error is due to the uncertainty in 
the measured moments and the second error is due to 
uncertainty in the theoretical calculations Q. In the 
shape function scheme, using the exponential shape func- 
tion form 6|, m^sF) = 4.43t"o.o8 GeV/c 2 and = 
OA4t a ° 6 7 GeV 2 , with a correlation of -0.63. If the Gaus- 
sian shape function form were used, m^gF^ a- n( A Hn(SF) 
would increase by 0.13 GeV/c 2 and 0.01 GeV 2 , respec- 
tively. The spectra with the fitted parameters are com- 
pared to data in figure These results (without theory 
error) are then used to extrapolate the measured par- 
tial branching fraction from E 1 > 1.9 GeV to 1.6 GeV 
to allow comparisons to theoretical predictions. In the 
kinetic scheme B(B -> X S ^,E 1 > 1.6 GeV) = (3.94 ± 
0.31 ± 0.36 ± 0.21) x 10~ 4 and in the shape function 
scheme B(B -> X sl ,E 1 > 1.6 GeV) = (4.79 ± 0.38 ± 
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OAitolr) x where the errors are statistical, sys- 

tematic and model-dependence. The model-dependence 
is derived from the la error ellipse for the m^-p^ fit. The 
central value in the shape function scheme is reduced to 
4.55 x 10~ 4 if the Gaussian form is used. 

Finally the sample is divided into b and b decays us- 
ing the charge of the lepton tag to measure Acp{B — > 

X s+ dl) = n++n- T^Z wnere are the positively 

(negatively) tagged signal yields and 1/(1 — 2uj) is the di- 
lution factor due to the mistag fraction lj. A requirement 
2.2 < E* < 2.7 GeV maximizes the statistical precision 
of the measurement as determined from simulated data. 
The yields are N+ = 349 ± 48 and iV - = 409 ± 45. The 
bias on Aqp due to any charge asymmetry in the detector 
or BB background is measured to be — 0.005±0.013 using 
control samples of e + e _ — > AT 7 and B — > Xtt°,t]. The 
mistag fraction due to mixing is 9.3 ±0.2% [lj|. An addi- 
tional 2.6 ±0.3% mistag fraction arises from leptons from 
D decay tt^ faking /i , 7 conversions, 7r° Dalitz decay 
and charmonium decay. After correcting for charge bias 
and dilution A CP = -0.110 ± 0.115(stat.) ± 0.017(syst.), 
including multiplicative systematic uncertainties from 
the BB background subtraction (5.4%) and the dilution 
factor (1.0%). The model-dependence uncertainty due 
to differences in the B — > X s j and B — > X^j spectra is 
estimated to be negligible. 

In conclusion, the branching fraction and the energy 
moments of the photon spectrum in B — > X s j are mea- 
sured for E 1 > 1.9 GeV. The moments are consistent 
with previous measurements 0, 0, 0] and are used to 
extract values of mj and /z^ which are consistent with 
those extracted from semi-leptonic B decays j^. These 
measurements have been used to reduce the systematic 
error in the estimation of |14&| an d |Kii>| 0- The mea- 
sured branching fractions are in agreement with the SM 
expectation and previous measurements. The measured 
Aqp is also consistent with the SM expectation. 
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Appendix 

The correlation matrices for the statistical, systematic 
and model-dependence errors of the first and second 
moments are given in tables |n] IIHI and IIVI respectively. 
The matrices are symmetric so only the upper half is 
tabulated. The moments are measured for four energy 
ranges, 1.9,2.0,2.1,2.2 < E 1 < 2.7 GeV. 
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TABLE II: The correlation matrix for the statistical errors between the the measured moments. FM and SM denote first and 
second moment respectively. 







FM 


tti r ni it fih it 

FM FM FM 


SM 


SM 


SM 


SM 


Quantity Min. _E 7 


1.9 


2.0 2.1 2.2 


1.9 


2.0 


2.1 


2.2 


FM 


1.9 GeV 


1.0000 


0.5172 0.3548 0.2265 


-0.6110 


0.0077 


0.0821 


0.1008 


FM 


2.0 GeV 




1.0000 0.6838 0.4326 


0.2285 


-0.0008 0.1375 


0.1645 


FM 


2.1 GeV 




1.0000 0.6260 


0.4220 


0.5660 


0.1650 


0.1884 


FM 


2.2 GeV 




1.0000 


0.4528 


0.7568 


0.7383 


0.2113 


SM 


1.9 GeV 






1.0000 


0.4626 


0.3486 


0.2520 


SM 


2.0 GeV 








1.0000 


0.6966 


0.4887 


SM 


2.1 GeV 










1.0000 


0.6709 


SM 


2.2 GeV 












1.0000 



TABLE III: The correlation matrix for the systematic errors between the the measured moments. FM and SM denote first and 
second moment respectively. 







FM 


FM 


FM 


FM 


SM 


SM 


SM 


SM 


Quantity Min. _E 7 


1.9 


2.0 


2.1 


2.2 


1.9 


2.0 


2.1 


2.2 


FM 


1.9 GeV 


1.0000 


0.7841 


0.6076 


0.4622 


-0.7552 


-0.3862 


-0.1723 


-0.0942 


FM 


2.0 GeV 




1.0000 


0.8007 0.6311 


-0.2515 


-0.4753 


-0.2298 


-0.1451 


FM 


2.1 GeV 






1.0000 


0.7886 


-0.0320 


0.0202 


-0.2903 


-0.2115 


FM 


2.2 GeV 








1.0000 


0.0489 


0.1850 


0.2061 


-0.1001 


SM 


1.9 GeV 










1.0000 


0.4473 


0.2219 


0.1436 


SM 


2.0 GeV 












1.0000 


0.5003 


0.3342 


SM 


2.1 GeV 














1.0000 


0.6592 


SM 


2.2 GeV 
















1.0000 



TABLE IV: The correlation matrix for the model-dependence errors between the the measured moments. FM and SM denote 
first and second moment respectively. 







FM 


FM 


FM 


FM 


SM 


SM 


SM 


SM 


Quantity Min. S 7 


1.9 


2.0 


2.1 


2.2 


1.9 


2.0 


2.1 


2.2 


FM 


1.9 GeV 


1.0000 


0.9267 


0.9486 


0.8252 


-0.9057 


-0.9223 


-0.9234 


-0.7983 


FM 


2.0 GeV 




1.0000 


0.9576 0.8669 


-0.7452 


-0.8183 


-0.8087 


-0.8960 


FM 


2.1 GeV 






1.0000 


0.9540 


-0.7390 


-0.7889 


-0.7843 


-0.7752 


FM 


2.2 GeV 








1.0000 


-0.5378 


-0.5857 


-0.5788 


-0.6023 


SM 


1.9 GeV 










1.0000 


0.9650 


0.9824 


0.6983 


SM 


2.0 GeV 












1.0000 


0.9810 


0.8035 


SM 


2.1 GeV 














1.0000 


0.8057 


SM 


2.2 GeV 
















1.0000 



